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Abstract 

We determine the patterns of spontaneous symmetry breaking in strong-coupling lattice QCD in 
' a fixed background baryon density. We employ a next-nearest-neighbor fermion formulation that 

possesses the SU(Nf) x SU (Nf) chiral symmetry of the continuum theory. We find that the global 
symmetry of the ground state varies with Nf and with the background baryon density. In all cases 
the condensate breaks the discrete rotational symmetry of the lattice as well as part of the chiral 
g | symmetry group. 

(N 

CM 
> 

CO 
O 

o 

CO 

o 

<D 
> 

X 
S3 



1 



I. INTRODUCTION 

The study of quantum chromodynamics at high baryon density, nearly as old ^ as 
the theory itself, has gained impetus in recent years with new interest in the idea of color 
superconductivity (CSC) [1, H HI HI ■ Recent work has led to a rich phase structure, including 
the possibilities of color-flavor locking |(| and crystalline superconductivity 0]. For a review 
see |8(. 

All this work depends on effective theories derived from (or at least motivated by) weak- 
coupling QCD. The running coupling, however, will become weak only at high densities; 
in fact it turns out that reliable calculations demand extremely high densities 0, llpj . If 
any of the predictions for moderate densities are to be believed, they must be confirmed by 
non-perturbative methods, or at least by models that incorporate QCD's strong-coupling 
features. 

In a previous paper jll| we constructed a framework for introducing baryons to lattice 
QCD at strong coupling. Though the lattice theory is a way of defining QCD at the most 
fundamental level, this process requires weak coupling; the strong-coupling theory may be 
regarded as an effective theory at large distances. It displays the non-perturbative effects 
of color confinement and spontaneous breakdown of chiral symmetry (for suitably chosen 
fermion formulations). Our framework is based on the Hamiltonian approach. We use strong 
coupling perturbation theory to write an effective Hamiltonian for color singlet objects 
[l2l Fill . At lowest order we get an antiferromagnetic Hamiltonian that describes meson 
physics with a fixed baryon background distribution. (Baryons move only at higher order.) 
We study this Hamiltonian through its path integral, which takes the form of a nonlinear a 
model. 

The strong-coupling theory has no color degrees of freedom. Its properties can be deter- 
mined through study of its global symmetries. The global symmetry group of the action de- 
pends on the formulation of the lattice fermions. For naive, nearest-neighbor (NN) fermions 
the symmetry is U(N) with 

N = 4N f , (1) 

realized on the quark spinors by combining the Dirac index with the flavor index. This 
too-large symmetry is indicative of the doubling problem of naive fermions 14]|. Adding 



longer-range terms to the fermion kernel can reduce this artificial symmetry. We add a next- 
nearest-nearest neighbor (NNN) term inspired by the SLAC fermion formulation (12I . 
In weak coupling, of course, this theory is still doubled; in strong coupling, there are no 
apparent ill effects of this doubling, and the U(N) symmetry is broken to the U(Nf) x U (Nf) 
symmetry of continuum QCD. 1 



II. EFFECTIVE HAMILTONIAN AND a MODEL 

In we studied the NN theory only, and showed that its global U(N) symmetry is 
spontaneously broken to a subgroup that depends on the baryon number. In this paper we 



The axial U(l) symmetry of the lattice theory is of course not present in continuum QCD, but it is 
inevitable on the lattice if we start with a local, chirally symmetric fermion theory 14]. We can mend 
our effective theory by hand, by adding new terms derived from an 't Hooft instanton vertex. 
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add the NNN terms to the Hamiltonian and thus study a theory with the same symmetry 
as the continuum theory. The effective Hamiltonian in strong coupling is 11 



H = JiT, QIQIh + J 2 E ««• (2) 

ni,r? 711,77 

Here Q\ are U(N) charges at site n, with r\ = 1, ... , N 2 . This Hamiltonian moves mesonic 
excitations around the lattice, leaving the baryon density fixed. The states at n comprise a 
U(N) representation whose Young tableau has N c columns. The number of rows to depends 
on the baryon number B at n according to 

to = B + 2N f . (3) 

The sign factors are given by 

sjf = 2TiM v a i M v a i . (4) 

Here M v are the matrices of the fundamental representation of the U (N) algebra, normalized 
in the usual way, TtM v M v ' = \8 m >. The matrices oti are the 4x4 Dirac matrices times the 
unit matrix in flavor space. 

The NN term in the Hamiltonian is that of a U(N) antiferromagnet, while the NNN terms 
break the symmetry to U(Nf)i, x U(Nf)n. If one derives the fermion Hamiltonian by trun- 
cating the SLAC Hamiltonian, then both couplings J\ and J 2 are positive, and J 2 = Ji/8. 
If we argue, however, that the strong-coupling Hamiltonian is derived by block-spin trans- 
formations applied to a short- distance Hamiltonian, then we cannot say much about the 
couplings that appear in it. We will assume that couplings in the effective Hamiltonian fall 
off strongly with distance, that is, < Ji <C J\. 

We use spin-coherent states [l6| to write the partition function for the effective Hamil- 
tonian (J2J). This is the path integral for a Euclidean nonlinear a model. The a field at 
site n is an iV x iV hermitian, unitary matrix that represents an element of the coset space 
U(N)/[U(m) x U(N — m)]. It can be written as a unitary rotation of the reference matrix A, 

a n = U n AUl, (5) 

where 

A =(V-iL)- (6) 

and U n G U(N). 

The action of the a model is 



S = ^ / dr 



- J2 Tr MJ ] n d T U n + ^ Yj Tr (VnVn+i) + E Tr {^n^iOn^lOli) 



(7) 



The NN term is invariant under the global U(N) transformation U — > VU, or a —> VaV^ . 
The NNN term is only invariant if V^a^V = oti for all i. This restricts V to the form 

V = exp [i (6 a v + 75 ^) X a ] , (8) 

where A a are flavor generators. This is a chiral transformation in U(Nf) x U(Nf). 
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The NNN term couples (discrete) rotational symmetry to the internal symmetry, viz. 

a n — > R^cTniR, ri = Tin. (9) 
Here 72. is a 90° lattice rotation and R represents it according to 



R = exp 







)1 




V° 


a i J. 



® ljv r 



(10) 



III. NEAREST-NEIGHBOR THEORY 



The overall iV c factor in Eq. (JJJ) allows a systematic treatment in orders of 1/N C . In 
leading order, the ground state is found by minimizing the action, which gives field configu- 
rations that are r independent and that minimize the interaction terms. For the NN theory, 
minimizing the single link interaction 



E = y Tr OlOa 



(11) 



allows us to construct the vacuum by placing a\ and cr 2 on the even and odd sites. 

We impose a uniform baryon density, B n = B > 0, on the effective Hamiltonian by 
setting a fixed value of m > 2Nf on every site. To minimize the single-link energy (fTTf we 
first choose a basis where 

a 1 = A=( 1 ™ , V (12) 



N-m 



The analysis in [11| then shows that 02 can take any value of the form 



C"2 




(13) 



where the m x m submatrix can be chosen freely in the submanifold 

U(m)/[U(2m — N) x U(N — m)] according to 



a 



(m) _ jj(m) ^(m)jj{m)] ^ 



with 



^-2m-N 
ljV-m 



(14) 



(15) 



and £/( m ) g U(m). As mentioned, we construct a ground state of the infinite lattice by 
replicating o\ and 02 on the even and odd sites of the lattice. Thus while all the a n on the 
even sites point to A, on the odd sites each a n wanders independently in the submanifold. 
This classical ground state has a huge degeneracy, exponential in the volume. 

In Ref. |ljj we showed that the 0(1/N C ) fluctuations generate a ferromagnetic inter- 
action am ong the odd sites, causing them to align to a common value ("order from 
disorder" [17[). This is a Neel structure, with two sublattices. The even sites break 
U(N) to U (m) x U (N — m) and then the odd sites break the symmetry further to 
U(2m — N) x U(N — m) x U(N — m). 
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IV. NEXT-NEAREST-NEIGHBOR THEORY 



Now we add the NNN interactions to the effective action. At the cfassical levef, they do 
not by themselves remove the classical degeneracy of the NN theory. We have to introduce 
the 0(1/N C ) fluctuations first in order to stabilize the Neel ground state of the NN Hamil- 
tonian. Hence we assume 1/N C > J2/J1, and treat the NNN interactions as a perturbation 
that lifts part of the (global) degeneracy of the 0(1/N C ) ground state. 

We begin, then, by assuming a Neel ansatz that minimizes the NN term in the action 
(J7J). The NNN term acts within each of the two sublattices. Writing a e>0 for the sublattice 
fields, the NNN contribution to the energy per 2x2x2 lattice cell is 

E nnn = y \^aOLiO a OL^ . (16) 

a=e,o i 

Here a a is a global unitary rotation of the solution to the NN theory given by Eqs. fT2|) - (fH 

We can find a lower bound for E nnn . Writing in each term = a a and £!f = a a i<J a ot\^ 
we have 

E nnn = ~ Tr E°Ejf. (17) 



Note that S 12 are unitary rotations of the reference matrix A. Each term in Eq. ()17j) may 
be bounded from below by allowing these unitary rotations to vary independently over the 
entire U(N) group. This is just the minimization problem posed in Eq. (fTT|) above. The 
solution is given by Eqs. (|T2^) - (fT3|) . whence the bound 

£„„„ > 3J 2 (4m-3A). (18) 

We minimize E nnn by writing an ansatz for a a that saturates the lower bound. At this 
point we choose to work in a basis where 75 is diagonal, 

- - ( Y -L ) • (l9) 

Our ansatz is 






UA a lP = U -l N -m I U ] . (20) 



1 



N—ra 

(Note that A Q is a rotation of A e .) This is a global rotation (via U) of a particular configu- 
ration that minimizes the NN energy, as seen above. We further suppose that U takes the 
form 

(21) 

with u e U(N/2). 

The minimization of the energy via the ansatz is a problem of vacuum alignment 18]. We 
begin with the unperturbed NN problem, where the symmetry group G = U (N) is broken to 
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H = U(2m — N) x U(N — m) x U(N — m). Our reference vacuum is given by a ep = A 6)0 , 
giving a specific alignment of H as the invariance group of this vacuum. We determine the 
rotation matrix U [within the ansatz ([21)1] that minimizes the energy of the perturbation 
E nnn , which we write as 

E mm = ^ E Tr AoSiAaaJ, (22) 



where ot.i = U^aciU = a]. 

We denote the generators of if collectively as T and the remaining generators of G as X. 
The T matrices commute with both A e and A D , while the X matrices do not. The symmetry- 
breaking term in the energy is E nnn , given by Eq. (|22|) in terms of the rotated Hermitian 
matrices «j. We project each oti onto the T and X subspaces, giving the decomposition 



a; 



— T — X 
<X; +01; . 



(23) 



Using the invariance of A e o under T, 



(24) 



and the orthogonality of the T and X subspaces, 



Tr 



(25) 



we have 



En 



J 2 



2 f 



^Tr [(af) 2 + A a afA a «f 



(26) 



Following the block form of A e)0 , we divide the broken generators X into 3 sets, denoting 
them as X a with a = 1, 2, 3. Their structures are respectively 



/ X x 



x 2 



ij , 
\ 00/ 




(27) 



Writing af = J2a a i a ; t ne following relations can be proved easily: 



A^ctf 1 ] = {A^af 1 } = 
{Ai,«f 2 } 
{Ai,«f 3 } 








x 2 

i 



A 2 , a 
{A 2 ,«f 3 } 




0. 



Using these together with the further orthogonality conditions, 



(28) 



Tr 



a 7^ b, 



we bring Eq. (|26|) to the form 

E nnn = J 2 J2^ \(ajf - (af 3 )' 



(29) 



(30) 
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The rotation U given in Eq. (J2~T|) saturates the lower bound for the energy. We will 
proceed to prove this for the case m > 3N/4. In view of Eq. (}2~T]) we can write a, in the 
form 

*=(* 4 o)' (31) 

where <7j = wVjW. It is straightforward to check that for m > 3N/4 we have 

(32) 
(33) 

with a\ composed of the first (2m — 3N/2) columns of <7j, 

( -n for q = l,...,2m-3N/2, 

M«= else. (34) 




The energy is 



2m-37V/2 

^nnn = 2J 2 El(^) P9 | 2 = 2J 2 E E = 3J 2 (4m - 3iV), (35) 

i,pq i q=l 



which is exactly the lower bound. 

According to Eq. (J33J), the bound is saturated for any u G U(N/2). Different choices of u 
are not in general related by transformations of the U(Nf) x U(Nf) symmetry group. There 
is thus an accidental degeneracy of the vacuum when the NNN term is treated classically. 
This degeneracy is presumably lifted by fluctuations. 

The simplicity of this calculation depends on the assumption m > 3N/4. For m < 3N/4, 
both af A and af are nonzero. Moreover the index structure of the projections is more 
complex. 2 Therefore in these cases we resort to numerical minimization of Eq. (J22j) over 
u G U (N/2). In each case we find that the ansatz (j2(Jj) -(|21j ) yields a minimum that saturates 
the lower bound Again, there is the possibility of accidental degeneracy. 

Upon calculating the a fields using Eqs. (|2U|l - (j2T|l . it is straightforward to ascertain the 
symmetry of the vacuum. The U(Nf) x U(Nf) generators that commute with both a e 
and o form the unbroken subgroup of the NNN theory. The rest are broken generators 
that correspond to Goldstone bosons. We summarize our results in Table HI and we note the 
following: 

• In cases of accidental degeneracy, we show the largest unbroken symmetry attainable. 
For m > 3N/A (i.e., B > Nf), this comes of the choice u = In/2- For m < 3N/A our 
numerical work cannot rule out vacua with yet larger symmetry. 

• Since the baryon background is fixed, we cannot tell whether the U(l) corresponding 
to baryon number is broken. [The U{1)b group acts trivially on Q 1 ^ and on a n .) 

• For each value of Nf, the case B = 2Nf is a completely saturated lattice. Each site is 
in a singlet under the chiral group, and there is no spontaneous symmetry breaking. 



2 An exception is the m — N/2 case (B — 0), which was solved in 
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TABLE I: Breaking of SU(Nf)L X SU(Nj)r x U(1)a for all baryon densities (per site) accessible 
for Nr < 3. Results for B = are from 



N f 


\B\ 


Unbroken symmetry 


Broken charges 







— 


1 


1 


1 


— 


1 




2 


U(1)a 










SU(2) V 


4 




1 


U(l)i 3 


6 


2 


2 


SU(2) V 


4 




3 


u(i) h 


6 




4 


SU{2) L x SU(2) R x U(1) A 










SU(3) V 


9 




1 


U(1) Y x SU(2) V 


13 




2 


U(1) Y 


16 


3 


3 


SU(3) V 


9 




4 


U(l) h x U(V)y 


15 




5 


U(l) h x U(1) Y x U(\) A , 


14 




6 


SU(3) L x 5C/(3) R x U(1) A 






The axial U(l) is not a symmetry of the continuum, and must be broken by hand. 
Wherever it appears in Table H] it should be neglected, whether as an unbroken sym- 
metry or as a broken charge. 



The U(1)a' appearing in the table for Nf 
rather is generated by 75 £g> A' with 



3 is not the original U{1) A group but 



A' 



( I o\ 





(36) 



This is the only case where an axial symmetry survives spontaneous symmetry break- 
ing. If U{1)a is broken by hand, so is U(1) A >- 

For all nonzero densities short of saturation the vacuum breaks rotational invariance. 
This is easily checked for m > 3N/A by noting that the ansatz for <7 ej0 fails to commute 
with the rotation operator (|10p. For the other cases this is easily checked numerically. (In 
some cases a discrete symmetry around the z axis remains unbroken.) Since this is not 
a continuous symmetry it will not give rise to additional Goldstone bosons. The broken 
rotational invariance will of course affect the excitation spectrum. In particular, whereas 
the NN theory possesses excitations with linear and quadratic dispersion relations .17], the 
NNN theory will produce interesting admixtures with anisotropic dispersion relations, like 
those seen in We defer discussion of the excitations to a future publication. 
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V. DISCUSSION 



In comparing our results to those of continuum CSC calculations, one must keep in 
mind that we study systems with large, fixed, and discrete values of B, rather than with 
large, continuous fi. Moreover, we use large- N c approximations which necessarily ignore the 
discrete properties of the SU (N c ) group that are essential to baryons. Quantum effects at 
finite N c , treated correctly, should yield effects that are not accessible through the 1/N C 
expansion. 

The values of (N c , Nf) that are of interest for CSC are N c = 3 and Nf = 2 or 3. In the 
two-flavor case the favored qq condensate is a flavor singlet and a color triplet, so that while 
color is partially broken, chiral symmetry is unbroken. We do not see this for any density 
at Nf = 2. Plainly our results are due to a qq condensate; whether there is a qq condensate 
as well cannot be ascertained. 

For Nf > 2 the situation is similar. Schafer |2Q] has considered the color-flavor structure 
of the condensates that arise for N c = 3 and Nf > 3, and he has found that both color 
and flavor are partially broken, with a condensate that locks one or more subgroups of the 
flavor group to the color group. Since we work at large N c , we should stand the argument 
on its head. A plausible qq condensate would lock successive subgroups of the color group 
to the flavor group and hence to each other, leaving unbroken the diagonal SU(Nf)L + R + c 
and some leftover color symmetry. Judging by the global symmetry of the vacuum, perhaps 
we see this for (Nf — 3, B — 3). The other cases could conceivably arise from a combination 
of qq and qq condensates, but whether the latter actually occur is an open question. 

Finally we note that according to diagrammatic power-counting arguments 21, 2^1, CSC 



should disappear in the 't Hooft limit (N c — > oo, g 2 N c fixed). We do not strictly work in 
this limit, since the large- N c approximation is applied only to the effective strong-coupling 
Hamiltonian, after g 2 has disappeared into setting the energy scale. 
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